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AMBIENT  SEA  NOISE  DIRECTIONALITY: 
MEASUREMENT  AND  PROCESSING 


INTRODUCTION 


In  the  study  of  ambient  noise  in  the  ocean,  the  spatial  coherence 
of  the  noise  field  can  be  characterized  by  the  angular  distribution  and 
strengths  of  the  noise  arrivals.  Such  a description  of  the  noise  field 
is  generally  very  useful  in  analyzing  the  causes  of  ambient  sea  noise 
and  its  relation  to  the  ocean  environment.  Furthermore,  ambient -sea- 
noise  directionality  data  provide  the  basic  information  in  designing 
the  receiving  array  used  in  underwater  communication  systems. 

A considerable  amount  of  work  concerning  ambient-sea-noise  direc- 
tionality  has  been  conducted  since  the  early  part  of  1960.  As 

noise  directionality  is  not  a quantity  which  can  be  measured  directly, 
many  different  processing  methods  have  been  employed  to  extract  ambient- 
noise-directionality  data  from  the  output  of  an  array  of  receiving  hydro- 
phones. The  theoretical  analysis  of  the  relationship  between  the  array 
output  and  noise  field4 >5, 6 indicates  that  the  exact  noise  field  (unique 
solution  for  the  mathematical  problem)  is  difficult  to  determine,  and 
only  limited  information  about  the  noise  field  (principal  solution)  can 
be  recovered  by  an  array  with  a finite  aperture  size. 

Techniques  developed  to  estimate  the  noise  field  are  based  on  either 
the  array  beamformer  output^>3,7  or  cross-spectral  matrix  between 
its  hydrophone  outputs.®  Because  of  the  complications  in  deriving  the 
processing  algorithm,  discussions  of  those  techniques  are  generally  con- 
fined to  an  equally  spaced  hydrophone  line  array.  Utilization  of  a 
deconvolution  technique  for  nonuniformly  spaced  arrays  has  been  reported 
recently  in  the  literature.®  The  analytical  properties  of  such  a spe- 
cific algorithm  have  not  been  elaborated  on  yet. 

On  a separate  front,  research  efforts  have  been  directed  to  the 
improvement  of  array  design,  the  measurement  tool  for  ambient-sea-noise 
t directionality.  Various  approaches,  some  of  which  are  extensions  of 

| knowledge  accumulated  in  fields  such  as  radio  astronomy  and  geophysics, 

t have  been  devised  to  optimize  array  performance  so  that  the  array  output 

| can  be  related  directly  to  the  noise  field.  Those  approaches  usually 

l involve  computation  of  a set  of  shading  factors,  amplitude  and  phase, 

applied  to  array  hydrophones  to  maximize  certain  array  performance  param- 
eters: directivity  index, 10*11  main-lobe  to  side-lobe  ratio, 12* 13  or 

1 high-order  pressure  gradient. 


1 


Depending  on  which  array  parameter  has  been  optimized,  the  array 
beamformer  output  will  have  least  interference  from  its  side-lobe  res- 
ponse; hence,  its  measurement  can  be  more  closely  related  to  the  noise 
field  in  the  direction  in  which  the  array  is  steered.  Since  the  noise 
field  is  not  known  beforehand,  the  optimization  of  array  performance 
has  to  be  based  on  an  assumed  noise  field,  generally  isotropic,  which 
is  not  necessarily  true  in  reality.  Such  a shortcoming  is  eliminated 
only  after  the  optimum-adaptive  beamforming  technique^  is  adopted. 
However,  the  processing  system  for  carrying  out  this  technique  is  com- 
plicated.^ in  addition,  even  if  the  array  gain  is  optimized  with  this 
processing  method,  the  width  of  the  array's  main  lobe  varies  with  steer- 
ing direction  and  depends  on  the  noise  field  it  intends  to  measure. 

Thus,  further  interpretation  or  processing  is  needed  to  express  the  meas- 
urement results  in  units  of  noise  level  per  steradian. 

| 

This  report  presents  the  results  of  an  investigation  of  measurement 
and  processing  methods  for  ambient-sea-noise  directionality  used  in  a 
receiving  array.  The  purpose  of  such  a study  is  to  emphasize  the  con- 
nections among  the  various  processing  techniques  and  to  establish  a 
common  reference  for  noise-directionality  data.  In  order  to  reconcile 
the  theoretical  analysis  and  experimental  measurement  of  ambient-sea- 
noise  directionality,  the  basic  relationships  between  array  parameters 
and  the  noise  field  are  reexamined.  Discussions  are  directed  also  to 
the  practical  problems  occurring  in  the  real  ocean  environment,  such  as 
stationarity  and  homogeneity  of  the  noise  field. 

Derivation  of  ambient-sea-noise  directionality  in  terms  of  the  meas- 
urable quantities  of  the  receiving  array  are  followed  through  in  great 
detail  in  this  report.  In  this  way,  a proper  physical  interpretation 
j can  be  applied  to  the  final  processed  ambient-sea-noise  directionality 

data.  the  limitation  of  the  numerical  computation  in  carrying  out  data 
reduction  has  been  addressed  with  consideration  for  availability  of 
, existing  computer  algorithms.  Although  the  overall  analysis  presented 

in  this  report  can  be  directed  to  any  type  of  array  used  for  ambient-sea- 
noise  directionality  measurement,  only  a vertical  AUTOBUOY  directional 
arrayl?  is  applied  as  an  example  in  the  computation  of  data  reduction. 
This  is  due  mainly  to  the  availability  of  field  data  for  the  AUTOBUOY, 

< so  that  comparisons  can  be  made  between  a simulated  noise  field  and  a 

real  noise  field  obtained  from  an  actual  ocean  environment.  It  is  hoped 
that  this  report  can  bridge  some  of  the  gaps  in  the  interpretation  of 
measurement  results  and  the  understanding  of  the  noise  field. 
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NOISE  FIELD  AND  ARRAY  PARAMETERS 


FUNDAMENTAL  RELATIONSHIPS 


Spatial  Correlation  of  Noise  Field 


A plane  wave  of  amplitude  n and  angular  frequency  u can  be  expressed 
mathematically  as 


ne 


-j(ut  - k*r) 


(1) 


where  k is  the  propagation  vector,  with  a magnitude  of  2r/X,  and  r is  a 
direction  vector  with  reference  to  a particular  set  of  coordinates. 


sponse 


A receiving  hydrophone  placed  at  s^  from  the  origin  will  have  re- 


= ne 


jk*s( 


(2) 


where  unit  gain  has  been  assumed  for  the  hydrophone  sensitivity  and  the 
time-dependent  factor  e‘Jwt  is  dropped  for  simplicity.  The  expression 
for  e in  equation  (2)  can  be  considered  as  a vector  (rotating  vector 
or  phasor)  whose  magnitude,  n,  and  phase  angle,  k*s£,  contain  the  infor- 
mation on  the  propagation  direction  of  this  plane  wave. 


Since  the  reference  of  the  phase  angle  is  chosen  arbitrarily,  meas- 
urements of  this  plane  wave  at  two  locations  at  least  are  needed  to 
retain  the  phase  angle  information.  As  shown  in  figure  1,  for  two  hydro- 
phones located  at  s ^ and  sm,  the  phase  angle  information  can  be  recovered 
simply  by  taking  the  cross  product  of  the  hydrophone  outputs: 

e e * = n2e^k  Sm^ 

l m 


= n2ejk'd*m 


(3) 


where  * denotes  complex  conjugation. 


When  the  noise  field  is  the  result  of  signals  arriving  from  various 
angular  directions,  both  n and  k will  be  functions  of  angle  n,  and  the 
total  response  for  the  cross-product  of  hydrophone  outputs  will  be 


e„e  * 
l m 


= /n  aca)ejtCn),a*m, 


dft 


(4) 
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Figure  1.  Measurement  of  Noise  Field 


Since  equation  (4)  concerns  only  a single  frequency  ui  (narrowband 
signal),  it  can  be  written  more  explicitly  as 


YfuO^m  = / NC^e^.^d,  , 


(5) 


where  y(w)£m  is  the  cross  spectrum  of  hydrophone  outputs  a and  m separa- 
ted by  d^m,  and  N(fi,io)  is  the  power  spectrum  of  the  noise  field  arriving 
from  direction  £1. 


For  noise  directionality  measurements,  an  array  of  several  hydro- 
phones will  be  employed.  Thus  different  values  of  will  be  obtained, 
depending  upon  the  hydrophone  spacing  in  the  array.  The  cross  spectrum 
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Y£m  can  be  related  to  the  noise  field  as  expressed  in  equation  (5) . 
Nuttall^  has  pointed  out  that  the  noise  field  can  be  estimated  directly 
from  the  cross-spectral  matrix  } without  further  array  processing, 
and  he  gave  a detailed  discussion.  Derivation  of  those  relationships 
will  be  examined  again  here,  but  with  a different  approach.  This  is  not 
only  to  verify  the  previous  result  but  also  to  establish  a new  concept 
for  array  processing. 

For  an  array  of  M hydrophones,  immersed  in  a stationary,  homogeneous 

noise  field,  there  will  be  a maximum  of  L = 1 + independent 

measurements  of  y , if  only  a single  frequency  is  considered  (the  inde- 
pendent measurements  will  be  reduced  to  M for  an  equally-spaced  hydro- 
phone array).  Consequently,  there  will  be  L sets  of  simultaneous 
equations  of  the  type  shown  in  equation  (5) . Using  matrix  notation 
(labeled  with  underline,  number  of  lines  will  designate  the  dimension  of 
the  matrix),  we  find 


71,1 

e^-dl.l 

71,2 

e j k * d x >2 

D = 

Y£,m 

Y(M-1),M 

ejk‘d(M-l),M 

Equation  (5)  can  be  expressed  as: 

R = / NDdfi  . (7) 


Recalling  that  jl  can  be  obtained  from  the  cross  products  of  array 
outputs  and  that  £ can  be  determined  according  to  the  array's  structure, 
only  N is  unknown  in  equation  (7).  Basically,  an  integral  equation  in 
the  form  of  equation  (7)  (Fredholm19  equation  of  the  first  kind)  will 
have  an  infinite  set  of  solutions  for  unknown  N.  However,  because 
there  are  available  a finite  set  of  independent  measurements  with  the 
array,  a unique  solution  for  N can  be  deduced  from  information  about  R 
and  D.  The  transformation  method  will  be  used  to  derive  the  following 
relations . 


Consider  D to  be  a function  whose  bases  are  e-1  ^lm.  Apply 
Schmidt's  orthonormalization  process^O  to  D and  express  JD  in  new  bases 
as 


5 


T*c  .1 


where  the  elements  of  Q are  determined  by 


/u.  u . *dSl 

i 1 


0 for  i i j 


or. 


df! 


I 

— > 


= 1 for  i = j 


(9) 


and  the  double  subscript  Jim  is  replaced  by  k with  the  relation: 

k = l + m - 1 . (10) 

I 

Since  the  transformation  matrix  Q is  a constant  matrix  (all  its 
elements  are  constants  determined  in  equation  (9)),  it  can  operate  on 
both  sides  of  equation  (7);  the  result  is 

( 


i 

t 
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where  t designates  an  adjoint  operation  (changing  columns  to  rows  and 
conjugating  the  elements)  and  hj  is  a constant  resulting  from  the  opera- 
tion D+Q+Q.  Derivations  of  equations  (6)  to  (12),  which  are  just  the 
application  of  standard  methods  described  in  most  text  books, 19, 20  are 
carried  out  in  great  detail  here  for  the  purpose  of  clarifying  the  nota- 
tions used  and  injecting  some  physical  interpretation  into  the  problem. 


Transformation  of  the  base  functions  in  equation  (8)  can  be  visu- 
alized as  a way  to  express  D and  R in  different  coordinates  in  a multi- 
dimensional space  (Hilbert  space).  The  term  "space  vector"  is  purposely 
avoided  for  describing  the  function  R (underlined  in  present  notation) 
in  order  that  it  will  not  be  confused  with  the  physical  space  vector, 
such  as  s and  k (overbar  in  present  notation)  defined  previously. 

The  result  given  by  equation  (12)  can  help  to  give  some  understand- 
ing about  the  measurement  of  noise  directionality  with  an  array.  If  the 
element  spacing  of  an  array  were  chosen  such  that  h^  is  directly  related 
to  the  cross  product  of  the  hydrophone  output  y ^ , that  is. 


h. 

i 

N 


(13) 


then  the  noise  field  could  be  deduced  immediately. 


An  example  of  such  a case  is  a linear  array  with  its  hydrophones 
spaced  at  half-wavelength  intervals  (or  multiples  of  a half-wavelength). 
Even  if  the  element  spacings  do  not  form  an  orthogonal  basis,  the  noise 
field,  according  to  equation  (12),  is  still  a linear  combination  of  the 

original  bases,  (e  ^ ^i).  Since  h^'s  are  constants,  and  only  {k*d^} 
contains  angular  information,  the  angular  resolution  of  the  noise  field 
is  determined  by  the  array-element  spacing.  The  finest  angular  resolu- 


tion of  the  noise  field  expressed  by  equations  (12)  or  (13)  that  can  be 
achieved  will  be  related  to  e ^ where  |5l|  is  the  largest  array- 
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cos[kd  cos  (90  - a) ] = 


kdL  Sin  a = 3 

(14) 

,1  1 ,X  . 

Sin  “ = 3%>d£  = 6(dT} 

a = j(^)  (if  a < i)  . 

Similarly,  for  endfire  response  when  k is  parallel  to  5l,  the  maxi- 
mum magnitude  of  the  real  part  of  e dL  occurs  when 

cos [kdL (cos  6 - 1)]  Q = 0 = i . 

The  angular  deviation  6 for  a 3 dB  change  in  response  can  be  traced 


coslkd^l  - cos  8)]  = y 


kdL(l  - cos  6)  = y 


c°s  6 = 1 . I(A_) 
b - (if  e < h . 


Whether  a measurement  from  an  array  will  attain  such  an  angular 
resolution  for  a noise  field  depends  on  the  actual  distribution  of  the 
noise  field  and  the  uncorrelated  noise  in  the  measurement  environment. 
However,  the  resolution  criteria  suggested  by  Bracewel 1 and  Roberts, ^>4 

-—t—  in  the  broadside  direction  and  h—  in  the  endfire  direction,  are  not 
2dL  dL 

quite  correct.  They  apply  only  to  the  case  where  each  hydrophone  ele- 
ment has  equal  weighting  (shading)  factor.  Bases  other  than  those  shown 
in  equation  (13)  may  be  used  to  improve  the  angular  resolution  of  the 
processed  noise  field.  The  extension  and  discussions  will  be  presented 
later. 


Finally,  it  should  also  be  realized  from  equation  (12)  that  if  N is 
~ i k • d * 

formed  by  the  bases  (e  J x),  its  coefficients  h^  can  also  be  determined 
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5 ' I cifi  ' FA 


and  the  integrated  square  of  the  error  in  comparing  with  the  true  noise 


E = f (N  - N)2dfi 
L 

= f (N  - £ c.f.)2dfi 


By  setting  -= — = 0 to  minimize  the  error,  it  was  found  that 
dci 


f . = e 
1 


= e-jC‘ai 


C = A~  1 R 


Equations  (19)  and  (21)  then  agree  with  equation  (17).  According 
to  Bessel's  inequality, 20 


f ( N - ^ hie'^'ai)2dfi  >_  0 , 


/ (N  ' £+H)2dfi  >_ 0 . (22) 

The  integrated  squared  error  will  vanish  only  when  H is  a complete 
basis  set;  that  is,  a complete  basis  set  for  H is  chosen  to  represent 
N.  Elaboration  on  the  physical  interpretations  of  this  conclusion  will 
be  discussed  further  in  the  next  section. 

Array  Beamforming  Output 

The  beamforming  output  of  a receiving  array  has  been  used  as  a tra- 
ditional method  to  measure  the  noise  field  by  steering  the  main  lobe  of 
the  array,  either  mechanically  or  electrically,  in  various  directions. 
The  response  of  the  array  is  the  square  of  the  sum  of  each  hydrophone 


'*♦*  *•  TVW  * t<r  .* 
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output,  which  can  be  expressed  as: 


r M 
* /"  1 "l) 

it 


ej(k-s.  - ks.fik)  ^ > 


where  is  the  array  beamforming  output  for  the  array  steered  at 
(Wik)  are  the  shading  factors  applied  to  the  respective  array  elements. 
From  the  results  of  several  independent  measurements,  bj,  at  different 
steering  angles,  the  vector  form  of  equation  (23)  can  be  rewritten  as: 


where 


B = j NGdfi  , 


2,  W..ejCR's'i  - ksinl)  : 
i=l  n 

i?i  Wikej(K,ii  ' kSiV  ' 


^ W.Tej(t’5i  - ksiV  : 

_i=l  iJ 


Notice  that  equation  (24)  has  the  same  form  as  equation  (7).  Hence, 
the  same  method  used  to  solve  for  N in  equation  (7)  (discussed  in  the 
last  section)  can  also  be  applied  to  equation  (24)  to  deduce  N in  terms 
of  the  measured  beamformer  output  £ and  the  array  beam  pattern  £ at  the 
respective  steering  angles.  The  final  noise  field,  N,  will  be  expressed 
by  the  bases  of  £,  so  the  angular  resolution  of  the  computed  noise  field 
will  depend  on  the  beam  pattern  of  the  array.  Axelrod,  et  al.,3  derived 
the  same  result  by  a least-square  criteria  approach  and  described  in 
great  detail  its  application  to  real-field  data  processing. 


Examining  equation  (24A)  indicates  that  the  bases  of  £ are  a func- 
i k * s * 

tion  of  eJ  1 , which  is  related  to  the  array-element  spacings.  There- 
fore, further  simplification  of  the  equation  can  be  pursued  to  find  the 
direct  relationship  between  the  noise  field  and  the  array  element  spac- 
ings. Letting 
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equation  (23)  can  be  rewritten  as 

bk  - / N IwJeI" 

- / 


dn 


n (W^E  £+wk)  dn  . 


(25) 


(26) 


As  the  values  of  {W^}  are  determined  by  the  array  spacing,  steering 
angles,  and  shading  factor,  they  can  be  factored  outside  the  integral 
sign,  and  equation  (26)  becomes 


k£ 


bk  = 4+  / 

= \ f J NDdQV^ 


W.  +R  W. 
-k  = -k 


(27) 


■ 2 

i = l 


c.  y . 
l l 


where  the  definitions  of  FI,  I),  from  equation  (6)  and  (7)  are  used. 

Equation  (27)  gives  the  relationship  between  array  beamformer  output 
II  and  array  hydrophone  cross -spectral  matrix  R.  Since,  for  a given  array 

with  M hydrophones,  the  maximum  independent  set  of  R is  L = 1 + M_(M  ^ 1)  , 

the  maximum  number  of  independent  measurements  from  the  array  beamformer 
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output  for  the  same  array  will  certainly  not  exceed  L.  In  the  case 
where  the  number  of  measurements  obtained  from  the  array  beamformer  out- 
put equals  the  number  of  independent  Jl,  according  to  equations  (27), 

(7),  and  (12),  the  noise-  field  N deduced  on  the  same  bases,  (e  ^i), 
should  be  the  same  regardless  of  the  shading  factor  used  during  the 
beamforming  process. 

Such  a conclusion  can  be  deduced  easily  from  the  physics  of  the 
noise  field,  which  will  not  be  changed  just  because  certain  shading  fac- 
tors are  used  in  the  array  beamformer.  However,  when  the  number  of  meas- 
urements from  the  beamformer  outputs  is  less  than  L,  equation  (24)  must 
be  used  to  solve  for  the  noise  field.  Since,  under  this  condition,  not 
all  the  information  received  at  the  input  of  the  array  hydrophone  is 
utilized,  the  noise  field  deduced  via  equation  (24)  will  be  less  accu- 
rate in  general  than  that  given  by  equation  (7). 

REPRESENTATION  OF  NOISE  FIELD 

Coherent  Noise  Sources 


The  description  of  a noise  field  by  its  directionality  is  based  on 
the  concept  that  the  array  is  placed  in  an  infinite  space  (or  half-space) 
and  that  all  the  noise  sources  are  in  the  farfield  at  a distance  fur- 
ther away  than  L2/'A  (L  is  the  array-aperture  size,  A the  wavelength). 

Such  an  assumption  is  not  always  proper  for  ambient-sea-noise  measure- 
ments in  the  real-ocean  environment.  First,  propagation  of  an  acoustic 
wave  in  the  ocean,  particularly  over  long  distances,  is  influenced  by 
the  boundaries  of  sea  surface  and  ocean  bottom,  in  addition  to  the 
sound-veloci ty  profile  variation  along  the  propagation  path.  Even  for 
a single  noise  source,  the  wavefront  may  not  be  plane  and  many  paths 
may  be  taken  to  reach  the  measurement  site.  Second,  because  of  the  many 
interference  paths  taken  by  an  acoustic  wave  in  the  ocean,  the  noise 
field  (which  the  array  intends  to  measure)  is  not  only  nonstationary 
but  also  inhomogeneous.  Finally,  the  problem  with  nearfield  noise 
sources  is  that,  due  to  lower  propagation  loss,  they  are  strongly  coher- 
ent but  do  not  have  plane  wave  fronts.  Those  sources  will  not  be  des- 
cribed properly,  based  on  the  array  response  from  an  ambient-noisr- 
directionality  measurement.  Further  interpretation  and  representation 
of  a noise  field  in  this  type  of  situation  is  needed. 

The  conclusion  from  the  analysis  in  the  last  section  indicates  that 
the  basic  information  about  a noise  field  measured  with  an  array  is  con- 
tained in  the  input  to  the  hydrophone  and  the  correlation  of  the  hydro- 
phone responses.  Such  information  can  be  expressed  completely  for  an 
array  with  M hydrophones  as 
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(28) 


Here  the  omnidirectional  noise  levels  at  each  hydrophone,  diagonal 
terms  {-y^},  are  not  assumed  to  have  the  same  intensity.  Otherwise 
equation  (28)  can  be  reduced  to  equation  (7),  which  concerns  only  a homo- 
geneous noise  field.  Using  matrix  notation,  equation  (28)  can  be  rewrit- 
ten as : 


P = f N E E+ 


dfi 


(29) 


where  E^  has  already  been  defined  in  equation  (25). 


Expressing  the  measurement  results  in  the  form  of  equation  (25)  is 

based  mainly  on  choosing  the  particular  bases,  | Sm^ j , which  may 

not  be  the  best  representation  of  the  noise  field  that  the  array  is  meas- 
uring. To  examine  whether  further  relationships  between  which 

may  be  affected  by  the  multipath  of  noise  signals,  and  focusing  or  diffus- 
ing of  the  noise  field,  a mathematical  technique  of  diagonalizing  a 
matrix  is  employed  to  simplify  equation  (28).  Let 


A = S'1  P S = 


*1 

0 


0 


(30) 


= f N(S' 


N(S"lE)  (E+S)  dfi 


where  A is  a diagonalized  matrix  with  (Xm)  as  its  eigenvalues  and  S is  a 
matrix~composed  of  the  corresponding  eigenvectors. 
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The  physical  interpretation  of  this  operation  is  that  if  a new  meas- 
urement basis  (S'JE)  (E^S)  is  chosen,  the  intensity  of  principal  coher- 
ent noise  sources,  (XmT,  can  be  determined.  In  the  concept  of  a Hilbert 
space,  it  can  be  visualized  that  a set  of  coordinates  is  chosen  properly 
in  a multidimensional  space  so  that  the  space  vectors,  such  as  P,  can  be 
represented  simply.  According  to  equation  (30),  in  practice,  the  coher- 
ent noise  source  can  be  measured  directly  by  applying  its  corresponding 
shading  from  the  eigenvectors  of  S to  the  array  hydrophone-spacing  ma- 
trix E^.  It  should  be  realized  that  the  ability  of  an  array  to  measure 
such  independent  coherent  sources  depends  on  the  total  number  of  array 
elements,  M.  Although  the  derivation  presented  here  is  intended  to  ex- 
plore different  representations  of  the  noise  field,  this  technique  can 
serve  a very  practical  purpose  in  signal  detection  in  an  ocean 
environment. 

Slowly  Varying  Noise  Field 

The  angular  resolution  of  the  deduced  noise  field  from  array  meas- 
urements is  limited  by  the  array-aperture  size,  as  concluded  from  the 
previous  analysis.  This  is  because  all  the  available  information  is 
based  on  the  array-measurement  result.  However,  when  a noise  field  has 
low  variability  in  its  level  (very  weak  additive  independent  noise)  and 
there  is  little  disturbance  to  cause  uncertainty  in  the  measurement, 
the  noise  field  can  be  expressed  in  terms  of  other  bases  by  applying 
analytic  continuity  to  equation  (7).  For  example,  delta  functions  can 
be  employed  as  bases  to  represent  a noise  field  attributed  to  discrete 
sources : 


Substituting  equation  (31)  into  equation  (7), 


R = A • B 


(31) 


(32) 


where 
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Therefore  (b^)  can  be  solved  from  the  simultaneous  equation  in  (32), 
and  the  angular  resolution  of  the  noise  field  is  not  limited  by  the  ar- 
ray-aperture size.  However,  depending  on  the  narrowness  expected  in 
angular  resolution,  accurate  measurements  will  be  required  in  II,  which 
makes  equation  (32)  very  sensitive  to  numerical  error  (ill-conditioning 
of  equation  (32)).  This  will  impose  a limitation  on  the  practical  appli- 
cation of  equations  (31)  and  (32).  If  there  is  prior  knowledge  about 
the  noise  field  (unfortunately,  then  there  would  be  less  motivation  for 
making  a measurement),  perhaps  only  a certain  angular  section  of  the 
noise  field  would  require  fine  angular  resolution  for  examination.  This 
can  be  accomplished  easily  by  using  non-uniform  angular  resolution  in 
the  expansion  of  the  noise  field: 

N = £ b.Vffh  ± e.)  , (35) 

i 

where  V is  a window  function  defined  by  step  functions  u as 

V(ft.  ± e.)  = u(ft  -ft.  + c.)  - u(fi  - ft.  - e.)  . (36) 

111111 

Directionality  Ambiguity 

When  the  shortest  spacing  among  array  hydrophones,  ds,  is  larger 
than  a half-wavelength,  the  deduced  noise  field,  because  it  contains  a 

1i  k • d 

term  eJ  s according  to  equation  (12),  will  have  at  least  two  peaks  in 
the  angular  range  of  ft.  This  is  the  directional  ambiguity  caused  by 
using  such  types  of  arrays  in  ambient-noise-field  measurements.  Figure 
2 illustrates  this  problem.  The  cross  spectrum  (real  part)  of  two  hy- 
drophones plotted  versus  their  space  separation  has  a periodicity  of  X 
in  response  to  a wave  coming  from  the  endfire  direction.  For  an  array 

with  element  separations  larger  than  y,  the  sample  point  "a"  in  figure 
2 will  not  extract  all  the  information  (sample  rate  is  less  than  Nyquist 
. rate).  However,  if  the  general  directionality  of  the  noise  field  is 

known,  an  additional  data  point  can  be  interpolated,  such  as  "b"  or  "c" 

in  figure  2,  for  correlation  distances  less  than  y,  thus  avoiding  the 
directional  ambiguity. 

I 

» 
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PROCESSING  FOR  NOISE  DIRECTIONALITY 


DATA  COLLECTION 

The  primary  data  for  processing  for  the  noise  directionality  are 
either  the  cross-spectral  matrix  of  the  array  outputs  or  its  beamformer 
output,  as  has  been  discussed  in  the  previous  section.  However,  if  the 
stationarity  and  homogeneity  of  the  noise  field  has  to  be  analyzed 
beforehand,  the  original  information  is  contained  in  each  hydrophone  out- 
put, with  the  relative  phases  being  preserved.  This  requires  a simulta- 
neous recording  of  all  hydrophone  signals.  The  fluctuations  due  to 
array  motion  and  recording-tape  skewing  will  introduce  phase  and  ampli- 
tude errors  in  the  recorded  hydrophone  signal.  Methods  to  correct  these 
errors  should  be  devised  before  processing  the  data.  If  the  beamformer 
outputs  (real  positive  numbers)  or  the  cross-spectral  matrix  (complex 
numbers)  are  recorded  instead  of  the  hydrophone  signals,  information 
needed  to  verify  any  error  in  the  data  will  be  difficult  to  trace. 

By  examining  the  time  series  of  the  hydrophone  signals,  the  station- 
arity of  the  noise  field  can  be  analyzed,  so  that  proper  sampling  size 
and  rate  can  be  determined  for  processing  for  the  noise  directionality. 

As  the  array  occupies  a certain  physical  spacing,  the  problem  with  an 
inhomogeneous  noise  field  is  rather  bothersome.  Figure  3 illustrates 
an  example  where  measurements  obtained  for  the  cross  spectrum  not  only 
changed  with  hydrophone  separation,  but  also  varied  according  to  loca- 
tion. There  are  two  possible  ways  to  process  the  data:  pick  a particu- 
lar location  (line  CD)  and  read  all  correlation  values;  or  take  the 
average  value  across  all  positions  (line  AB) . Whatever  data  are  read, 
the  physical  interpretation  should  be  implemented  with  the  final  result 
of  noise  directionality  in  mind.  For  optimum  beamforming  processing^ 
in  estimating  the  noise  field,  the  inhomogeneous  noise  data  is  absorbed 
in  the  computation,  hence  the  physical  connection  is  also  lost. 


The  relationship  between  the  cross  spectrum  and  the  noise  field  is 
an  operation  of  linear  transformation.  Therefore,  the  results  of  sta- 
tistical studies  on  the  cross  spectrum  can  apply  also  to  noise  direc- 
tionality if  only  the  mean  value  is  desired. 

For  a large  densely  populated  array,  the  original  data  base  is 
large,  and  computation,  either  by  orthogonali zation  or  matrix  inversion, 
will  introduce  large  numerical  errors  in  addition  to  using  considerable 
computer  time.  A simplification  can  be  made,  if  the  noise  field  is 
nearly  stationary.  By  interpolating  data  points  at  locations  where  the 
base  functions  become  orthogonal,  the  noise  directionality  can  be  deduced 
from  those  interpolated  data.  If  the  smallest  separation  of  array 
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elements  is  larger  than  — , the  interpolating  process  may  eliminate  those 

noise  sources  whose  direction  is  close  to  the  endfire  direction.  Such 
problems  should  be  taken  into  consideration  when  the  original  data  are 
modified. 

If  the  noise  field  is  sure  to  oe  stationary  and  homogeneous,  or 
can  be  assumes  to  be  so,  simultaneous  measurement  of  the  cross -spectral 
matrix  is  not  necessary  for  deducing  the  noise  directionality.  Data 
from  towed  arrays,  or  a long-term  average  of  two  hydrophone  correlations 
at  various  separation  distances,  can  be  used  directly  for  processing  for 
noise  directionality.  Since  only  a sampling  of  cross-spectrum  data  is 
important  to  deduce  the  noise .directionality  in  a homogeneous  noise 
field,  an  equally-spaced  hydrophone  array  will  measure  redundant  data 
points.  An  analytic  approach  has  been  devised  by  Nuttall22  to  position 
the  array  elements  in  such  a way  that  equispatial  sampling  of  the  cross- 
spectral  matrix  can  be  achieved. 

The  angular  resolving  power  of  an  array  can  be  increased  if  the 
certainty  of  the  cross-spectrum  data  points  is  improved.  Hence,  infor- 
mation in  analyzing  data  fluctuation  and  their  statistics  will  be  useful 
in  final  processing  for  the  noise  directionality. 

COMPUTER  ALGORITHM 

Two  alternative  mathematical  procedures  can  be  used  to  process  the 
array  output  for  noise  directionality:  either  orthogonali zing  the  bases, 
or  solving  a set  of  simultaneous  equations.  There  is  a standard  numeri- 
cal method  for  orthonormalization  of  a set  of  functions2^  and  some 
available  program  routines  at  Naval  Underwater  Systems  Center  (NUSC)  can 
be  modified  easily  for  this  purpose. 2^  However,  this  method  involves 
numerical  integration,  where  accuracy  depends  on  the  data  segment.  Fine 
angular  resolution  in  noise  directionality  will  take  quite  a sizable 
computation  and  become  impractical. 

If  the  original  bases  of  the  array  are  not  orthogonal,  the  other 
approach  for  deducing  noise  directionality,  solving  a set  of  simultaneous 
equations,  will  be  easier,  since  the  matrix  element  in  A of  equation  (18) 
can  be  integrated  in  closed  form.  A subroutine  for  double-precision 
complex-matrix  operation  is  available  in  NUSC  Math  Pack.2^  ]n  solving 
simultaneous  equations,  the  major  operation  is  the  matrix  inversion, 

A 1 . The  accuracy  of  the  numerical  computation  depends  on  the  matrix 
"size  (related  to  number  of  array  hydrophones)  and  its  rank  (the  ability 
of  the  array  to  yield  independent  measurements,  related  to  element  spac- 
ing). The  error  introduced  during  the  numerical  computation  can  be 
estimated  from  the  eigenvalues  and  determinant  of  the  matrix. 26  When  the 
hydrophones  in  an  array  are  too  close  to  each  other,  the  matrix  A becomes 
ill-conditioned,  and  large  errors  will  exist  in  the  final  computation  of 
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noise  directionality.  To  improve  the  ill-conditioning  situation,  but 
degrade  the  resolution,  a small  fraction  of  deviation  can  be  added  in 
the  diagonal  terms  of  matrix  A. 

In  determining  the  principal  noise  sources,  resolution  of  the  eigen- 
values and  the  eigenvectors  for  matrix  P in  equation  (30)  is  needed.  A 
standard  computer  routine  can  be  found  Tn  the  NUSC  IMSL  Program  Library. 
The  information  contained  in  P can  also  be  used  in  conventional  beam- 
forming and  optimum  beamforming  processing. 21 


There  are  practical  problems  in  carrying  out  the  numerical  computa- 
tion for  finer  angular  resolution  of  the  noise  field,  as  discussed  ear- 
lier in  this  report.  First,  in  order  to  achieve  fine  resolution,  many 
samples  will  be  needed. from  the  cross- spectrum  function;  hence  the  size 
of  matrix  will  be  very  large.  Second,  because  of  fixed  array  aper- 
ture, sampling  points  may  be  too  close  to  each  other,  and  matrix  A can 
easily  become  ill-conditioned.  Meaningful  results  can  be  achieved  if 
there  are  only  a few  noise  sources  known  beforehand,  and  thus  dense 
sampling  will  not  be  necessary.  Therefore,  this  type  of  operation 
applies  only  to  some  particular  situations. 

Because  the  deduced  noise  field  is  based  on  the  summation  of  cer- 
tain base  functions,  if  those  base  functions  do  not  form  a sufficiently 
complete  set  to  express  the  noise  field,  the  final  results  of  the  sum 
will  not  be  exactly  equal  to  the  true  noise  field.  As  a result,  there 
will  be  some  oscillation  in  the  deduced  noise  field  data  (Gibbs  phenome- 
non) . 20  Smoothing  techniques,  weighting,  or  iteration  have  been  sugges- 
ted by  Nuttall8»22  for  the  processing  computations  to  reduce  such 
phenomenon.  However,  in  practice,  the  noise  field  is  not  known  before- 
hand. Therefore,  there  is  no  reference  to  the  method  for  applying  the 
smoothing  technique  to  improve  the  estimated  noise  field  (smoothing 
actually  degrades  the  angular  resolution  of  the  estimated  noise  field) . 


EXAMPLE 


| GENERAL 

The  analysis  results  in  the  previous  sections  apply  to  any  type  of 
array  in  general.  However,  the  examples  illustrated  here  have  been 
f limited  to  a 240  m vertical  array  operated  at  50  Hz.  This  is  the  array 

( developed  for  AUTOBUOY  as  a portable  system  for  ambient-sea-noise  direc- 

tionality measurements . 1^, 18  The  array  has  12  elements  spaced  unequally 
in  order  to  make  broad  frequency  band  (30  Hz  to  250  Hz)  measurements. 
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can  be  made  between  the  simulated  cases  and  the  real  ocean  environment 
in  processing  for  the  noise  directionality. 

Because  of  the  geometry  of  a vertical  array,  as  shown  in  figure  4, 
only  the  average  noise  in  the  cone  at  various  vertical  angles,  6,  can 
be  measured.  Hence,  the  general  relationship  of  the  noise  field  and  the 
array  parameters,  such  as  equation  (17),  can  be  reduced  to  a one-dimen- 
sional problem: 


/ 


R = / N (0) D sin  e d0  , 


0 


(37) 


where 


N (9) 


r2v 

■H- 


N (0 , 4>) dcj)  , 


(38) 


and  all  (k*d^)  can  simply  be  replaced  by  kdj  cos  6.  The  mathematical 
manipulation  can  be  made  ever 
to  u.  Equation  (37)  becomes 


manipulation  can  be  made  even  simpler  by  changing  the  variable  sin  (—  - 9) 


N(u)Ddu  , 


(39) 


and 


kd.  cos  9 — — kd. u . 
l i 


(40) 


NOISE  SIMULATION 

Figure  5 shows  the  estimated  noise  directionality  for  a broadside 
(9  = 90  deg)  noise  deduced  from  the  cross-spectral  matrix  of  12  hydro- 
phone outputs.  It  can  be  observed  that  there  are  small  oscillations 
around  zero  for  9 ^ 90  deg.  This  is  the  Gibbs  phenomenom  due  to  repre- 
senting the  noise  field  by  a finite  series  of  base  functions.  The 
cross -spectrum  matrix  P in  this  case  is: 


M 
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Figure  4.  Geometry  of  a Vertical  Line  Array 
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It  has  only  one  nonzero  eigenvalue,  12,  and  the  corresponding  eigen- 
vector is: 


S = — - 1 
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Figure  5a.  Cross  Spectrum 
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Figure  5b.  Estimated  Noise  Directionality 


Figure  5.  Broadside  Simulated  Noise  Field 
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This  set  of  values  can  be  used  as  shading  on  the  array;  then  the 
array  output  at  a steering  angle  of  90  degrees  will  be  equal  to  the 
noise  level.  Since  the  above  analysis  indicates  that  there  is  only  one 
noise  source,  we  can  assume 


N = a6  (0  - 60)  . 


According  to  equation  (33) 


4ikdx  cos  0 


r"  ikd 

i‘J  * 
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2 COS  0 


6(0  - 0O)  0 0 • • 0 sin  0 d0 


ikdio  cos  0 


^ikdj  cos  0O 


ikd2  cos  0( 


^ikdi2  cos  6o 


, ikd,  cos  0n  . „ 

1 = a e 1 u sm  0( 


sin  0O  • 


, ikd2  cos  60 

1 = a e ^ u sin  0O 


, ikdo  cos  0 n . . 

1 = a e 3 u sin  0O  . 


The  only  possible  solution  for  equation  (44)  is 


a = 1 


a - 1 
0 " 2 ; 
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hence 


N = 16(6  - ~)  . 

This  checks  with  the  original  assumption  of  a broadside  noise 
source . 

For  an  isotropic  noise  field,  the  cross  spectrum  can  be  expressed 
analytically  as 


Y . 


1 


sin(kdj) 

kdi 


+ jo  . 


(46) 


The  graphical  display  of  this  function  and  its  deduced  noise  direc- 
tionality are  shown  in  figure  6.  According  to  equation  (46),  the  cross 
spectrum  y^  equals  zero  at  multiple  half-wavelength  spacing.  As  the 
240  m array  corresponds  to  a length  of  15.56  half-wavelengths  at  50  Hz, 
the  eigenvalues  of  the  cross-spectral  matrix  in  this  case  will  have  at 
most  15  nonzero  values.  These  15  eigenvalues,  represented  as  15  inde- 
pendent noise  sources,  do  not  completely  describe  the  isotropic  noise 
field.  Therefore,  the  method  discussed  previously  is  not  applicable. 


The  cross-spectral  matrix  due  to  uncorrelated  local  noise  sources, 
such  as  system  and  turbulence  noise  generated  at  each  hydrophone  ele- 
ment, already  has  the  diagonalized  matrix  form.  Its  eigenvalues  will 
not  depend  on  the  length  of  an  array  and  the  locations  of  the  array 
elements.  For  an  array  shorter  than  a half-wavelength,  the  cross  spec- 
trum contributed  by  isotropic  noise  field,  sin(kd)/(kd)  will  be  approxi- 
mately one.  Thus  the  cross-spectral  matrix  (rank  one  matrix)  will  have 
only  one  dominant  eigenvalue.  However,  the  number  of  eigenvalues  due 
to  uncorrelated  local  noise  sources  will  not  be  changed.  Therefore,  in 
a super-directivity  array  application,  no  advantage  can  be  taken  by  ad- 
justing the  shading  of  array  elements  (changing  array  response  in  differ- 
ent-representation coordinates,  or  base  functions)  to  eliminate  the 
interference  caused  by  uncorrelated  local  noise  sources. 


The  effect  of  nearfield  noise  sources  on  the  array  response  is 
illustrated  by  a simulated  case,  shown  in  figure  7.  As  the  array  aper- 
ture is  240  m,  for  a 50  Hz  signal,  the  criterion  for  farfield  is 


a2 

f-  = 2020  m . 

A 

The  geometry  depicted  in  figure  7 gives  the  relationship  to  deter- 
mine the  cross  spectrum  among  the  hydrophones.  Figure  8 shows  the  esti- 
mated noise-field  directionality  based  on  the  cross  spectrum  computed 
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from  the  geometry  (range  = 120  m and  2020  m) . It  can  be  observed  from 
the  estimated  noise  directionality  that  the  array  has  no  way  to  detect 
a nearfield  source,  and  only  observes  a diffused  field.  Such  a short- 
coming can  be  eliminated  if  the  cross-spectral  matrix  is  diagonalized 
and  only  one  nonzero  eigenvalue,  12,  is  found.  The  corresponding  eigen- 
vectors (complex  number)  will  give  the  proper  shading  needed  for  array 
elements  to  focus  on  the  noise  radiated  by  this  nearfield  source.  If 
the  element  spacings  are  less  than  one  wavelength,  the  phase  angle  (ar- 
gument of  the  complex  number  shading)  can  be  used  to  reconstruct  the 
wavefront.  Hence,  by  geometry,  the  location  of  nearfield  noise  source 
can  be  estimated. 

SEA  DATA 

The  AUTOBUOY  array  has  12  unequally-spaced  elements,  since  the 
original  design  planned  to  cover  a broad  frequency  range. The  omni- 
directional noise  levels  obtained  from  the  12  hydrophones  for  various 
time  segments  are  plotted  in  figure  9.  These  curves  indicate  that  the 
noise  field  is  stationary  for  a sample  segment  greater  than  one  minute. 
However,  the  variation  in  the  noise  level  among  the  hydrophones  suggests 
that  the  noise  field  is  not  quite  homogeneous  (some  variation  may  also 
be  caused  by  system  noise,  such  as  that  caused  by  bad  connections). 
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The  noise  cross  spectrum  at  50  Hz  among  the  hydrophone  outputs  is 
displayed  in  figure  10  after  normalization.  Based  on  the  cross-spectral 
curves,  the  estimated  noise  directionality  is  plotted  in  figure  11.  The 
solid  line  in  that  figure  is  obtained  from  the  67  pairs  of  cross-spec- 
tral data.  However,  because  some  hydrophone  spacings  are  too  close, 

5 percent  deviation  has  been  added  in  the  diagonal  elements  of  matrix 
A for  avoiding  ill-conditioning.  The  dashed  line  is  based  on  the  inter- 
polated correlation  data  at  multiple  half-wave  spacings.  The  result 
obtained  from  the  optimum  beamforming  process^  is  included  in  figure 
12  for  comparison  (only  in  relative  scale).  It  can  be  observed  that  the 
general  intensity  distributions  are  about  the  same,  particularly  at  the 
peaks. 

Figure  13  depicts  the  computed  eigenvalues  from  the  noise  cross- 
spectral  matrix  measured  at  different  time  segments.  It  indicates  that 
there  are  at  least  9 major  dominant  noise  sources  within  a 10  dB  range. 
The  curve  falls  off  rapidly  from  #10  to  #12.  The  background  noise  is 
30  dB  down  from  the  major  noise  sources. 
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Figure  12.  Comparison  of  Optimum  Beamforming  and  Noise 
Field  Estimation  by  Transformation  Method 
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CONCLUSIONS 

From  an  analysis  of  the  noise  field  and  the  array  parameters,  it 
is  indicated  that  ambient-noise  directionality  can  be  derived  directly 
from  the  basic  information  contained  in  array  measurements,  the  cross 
spectrum  of  hydrophone  outputs.  If  the  same  array  is  used,  other  types 
of  measurements  should  also  reach  the  same  result  if  processing  methods 
are  correctly  applied. 

Theoretically,  the  noise  field  which  can  be  determined  from  the 
array  measurements  is  not  limited  by  the  aperture  size  of  the  array  if 
the  noise  field  is  composed  only  of  plane  waves.  However,  in  practical 
situations,  because  of  interference  due  to  uncorrelated  and  partially- 
correlated  noise  sources,  along  with  numerical  errors  introduced  during 
the  computation,  it  is  not  easy  to  achieve  such  a result.  The  most  re- 
liable solution  (principal  solution)  is  based  on  the  array  structure. 
For  a line  array,  computation  can  be  simplified  if  cross-spectral  data 
are  obtained  at  multiples  of  half-wave  spacings. 


The  basic  information  contained  in  the  cross  spectrum  of  an  array 
output  can  also  be  used  to  evaluate  the  major  coherent  noise  sources 
(eigenvalues).  The  corresponding  shading  factor  on  array  elements  (eigen- 
vectors) also  can  be  determined.  Hence,  the  optimum  response  of  an 
array  to  such  noise  can  be  achieved.  This  process  is  very  useful  to  de- 
tect nearfield  noise  sources,  but  cannot  be  implemented  in  the  noise- 
directionality  display. 

Since  the  cross  spectrum  of  the  array  output  gives  the  best  descrip- 
tion of  the  noise  field,  measurements  of  noise  field  should  be  designed 
to  sample  the  cross  spectrum  at  equally-spaced  intervals.  If  the  noise 
field  is  homogeneous,  then  the  array  spacings  need  not  necessarily  be 
equal.  Furthermore,  if  the  noise  field  is  stationary,  cross-spectral 
data  collected  at  different  time  intervals  also  can  be  used  for  process- 
ing for  ambient-noise  directionality. 


Deducing  the  ambient-noise  directional i 
data  is  a linear  transformation  operation, 
statistical  study  of  the  cross  spectrum  can 
directionality.  If  there  is  a high  degree  o 
data  (low  fluctuation),  fine  angular  resolut 
tionality  can  be  achieved  by  proper  use  of  p 
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